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ABSTRACT 

A degeneracy in strong lens model is shown analytically. The observed time delays and quasar image 
positions might not uniquely determine the concentration and the extent of the lens galaxy halo mass 
distribution. Simply hardwiring the Hubble constant {H^) and the cosmology (f2, A) to the standard 
ACDM cosmology values might not fully lift this degeneracy, which exists rigourously even with very 
accurate data. Equally good fits to the images could be found in lens mass models with either a mostly 
Keplerian or a flat rotation curve. This degeneracy in mass models makes the task of getting reliable 
Hq and A from strong lenses even more daunting. 

Subject headings: cosmological parameters — dark matter — distance scale — gravitational lensing 



1. INTRODUCTION 

One of the promises of gravitational lenses is to measure 
the Hubble constant Hq (Refsdal 1964) from the observed 
time delays among the images of a variable background 
quasar source lensed by a foreground galaxy. Given a 
model for the spatial distributions of the stars and dark 
matter in the lens galaxy, the time delay Atobs multiplied 
by the speed of light c is simply proportional to the ab- 
solute distances to the lens and the source, hence cAiobs 
scales with the size of the universe c/Hq. While the time 
delays are now routinely measured for many systems (see 
Schechter 2000), a reliable determination of ffp has been 
hampered to some extent by the intrinsic degeneracy in 
models of the dark matter potential of the lens (Williams 
& Saha 2000; Saha 2000; Zhao & Pronk 2001). The gen- 
eral trend is that a model with a dense dark matter halo 
gives a small Hq with 

TJoAtobs (X [2E„it(^/, zs) - ^.{Re) - Sh(i?£)] , (1) 
where Ecrit(2;,Zs) is the critical density for the lens at 
redshift zi and the source at Zs, and T,^{Re) or Sh(i?E) is 
the typical surface density of luminous or dark matter at 
within the Einstein ring Re (Falco, Gorenstein, & Shapiro 
1985; Kochanek 2002). 

Given that the value of Hq is now well constrained by 
other independent methods, such as the HST Key Project 
(Freedman et al. 2001), it is interesting to reverse the 
angle of the question, and use equation (1) to put more 
stringent constraint on the dark matter potential of the 
lens. More specifically in this Letter, we would like to ask 
the question: how narrow is the allowed parameter space 
for the lens dark halo which is consistent with a given set 
of images, time delays, cosmology and Hq? In the inter- 
est of clarity, we will consider only analytical lens models 
with a simplified geometry for a hypothetical image and 
lens system. We believe our arguments should apply qual- 
itatively to real galaxy lenses as well, and a more detailed 
apphcation to the quadruple system PG1115+080 is given 
in a follow-up paper (Zhao & Qin 2002). 

2. ANALYTICAL LENS MODELS FOR STARS AND HALO 
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For simplicity we assume the four images form a per- 
fectly symmetric Einstein cross with the time delay min- 
ima on the y-axis at ±Re (radian) from the lens center, 
and the saddle point images on the x-axis at zLqRe (ra- 
dian). For generality we rescale the images, taking the 
radius Re s^s unity in a new XY coordinate system, so 
the four images are given by 



-[0,±1], [A,Y],addic-[±g,0]. (2) 



This image configuration implies that a source {Xg, Yg) — 
(0, 0), exactly behind a spherical lens plus a linear external 
shear symmetric to the X and Y axes. 

All lensing properties can be derived from the time delay 
surfaces. For our model with spherical stellar lens (t>*{R) 
plus halo (j)h{R) and a linear shear of amplitude 71, the 
time delay contours are determined by a dimensionless 
time delay r(A, Y) given by 

(3) 

0,(i?)-0^(i?I,4) 



t{X,Y) ^t ■ uj{H^,^,zuZs)Re' 



li{X^-Y^) 



where R — VA^ + Y^ is the cylindrical radius, and 
uj{Hq, Q.,zi,Zs) is a constant containing all the dependence 
on the cosmological density parameter Q. and the lens and 
source redshifts z;, Zsi crudely speaking uj ^ Hq for typical 
lens and source redshifts in the ACDM cosmology. 

To illustrate the non-uniqueness, we need to find rea- 
sonable lens models with similar time delay surface. Let's 
consider the following example of lensing potential for the 
stars. 



0*(i?) ^ !^ln( 1 
a 



R° 



(5) 



where mo is the total stellar mass enclosed, a is the half 
mass radius, and 2 — a specifies the cuspiness of the stel- 
lar distribution. For the halo we use a nearly isothermal 
potential 

MR) = boR-- [(1, i?')min + HI,RX,.^ + ln\l,R^)^,^ 



(6) 
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where ^/bo is proportional to the terminal velocity of the 
halo rotation curve. 

Note that the function (j)h is smoothly connected at the 
two sides of the radius R = 1, which is chosen such that 
it is just outside the images. The parameter 5 is a dimen- 
sionless tunable parameter to adjust the halo contribution 
to the surface density at = 1. The surface density can 
be computed as 



k{R) = 1 



(7) 



so the densities for the stars and the halo are given by 



-h{R) = ^ 



2R [l,R^]r 



(8) 

Note that the stars have an inner cusp 2 — a, and the 
halo density is continuous across the break R= 1, and 
positive everywhere if bo/S > 2. 

By increasing S we can lower the mean densities at the 
images, hence creating the effect of a negative mass sheet. 
This will not affect the positions of the images, but can 
increase the time delay Arobs oc iloAiobs between the 
images. As we will see, this can result in a larger Hq 
to be consistent with the currently favored high value of 
Ho ~ 70kms"^Mpc"^ 

The deflection strength 



b{R) 



__ M,(< R) + Mh{< R) 



(9) 



dR R 

is effectively the rotation curve squared. A flat rotation 
curve corresponds to a constant deflection strength. For 
our model, the stellar and halo masses enclosed inside ra- 
dius R are given by 

moR° 



R^ + a"' 



Mh = boR-5 [(1, R^Uin + ln(l, R^)^^^] . 

(10) 

As we can see the stars converge to a finite mass mo at 
infinity and the halo dominates stars and approaches to a 
finite deflection (terminal velocity cx v^)- 



3. RESULTS 

Let's consider a typical lensing system in a standard 
ACDM cosmology with 

/fo = 70kms"^Mpc"\ O = 1 - A = 0.3 

with the lens and source redshifts 

zi = 0.5, Zs 

and a time delay Afobs 



2, 



(11) 



(12) 
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The observed time delay then set the following con- 
straint on the lens model, 

Arobs = t(0.5,0)-t(0,1) 



obs 



j-obs 



_ f 

saddle mimima 



-Ul 



(15) 
(16) 



(13) 

between the saddle point image at {qRE,0) and the min- 
ima image at {0,Re)- 

The above cosmology specifies the distances to the lens 
and the source, hence the time delay normalization con- 
stant 

Lo{Ho = 70,n = 1-A = 0.3, zi = 0.5, z^ = 2) ^ lOOkms^^Mpc 

(14) 



~ 28(1 - q^) day/D" x 100kms"^Mpc-^ ~ 0.22 (17) 

We solve for the parameters of the lens and the shear 
to reproduce the four images at {ztq,0), (0,±1) exactly. 
The images form at the extreme points of the time delay 
surface, hence we have the additional constraints 

— = 0, — = 0, at (X, Y) = (0.5, 0) and (0, 1). (18) 

It turns out that the allowed models follow a three- 
parameter scqueucx;. say (6, a, a). The parameters of a 
few representative models are listed in Table 1. The im- 
ages appear rigourously at the same locations for all four 
models, which are the minima and saddle points of the 
arrival time surfaces (cf. Figure la), and the relative time 
delay is identical as well for all four lens models (cf. Fig- 
ure lb). But the mass distributions in the four models are 
far from similar (cf. Figure 2a and Figure 2b): lensIII and 
lensIV correspond to systems with a finite stellar core with 
a very typical half-mass radius Re about half of the Ein- 
stein radius Re] lensIV is purely in stars, and has no halo. 
Both lensi and lensll correspond to systems with a strong 
stellar cusp dominating the halo at small radii; lensll has 
a bigger stellar component with a half-mass radius at one 
Re- 

3.1. The cause of degeneracy 

The above degeneracy is a variation of the well-known 
mass-sheet degeneracy. The latter implies that we can in- 
crease Hq by reducing the surface density k = ^ at the 

Einstein radius Re (cf. equation (1)), e.g., by either scal- 
ing down the stellar mass mo or scaling down the isother- 
mal dark halo mass b^r. But if we reduce the stellar mass 
mo and increase the halo mass b^r simultaneously such 
that we keep k^{Re) + Kh{RE) constant (between 0.10 to 
0.13, cf. Table 1), then /foAtobs will not change, only the 
terminal velocity of the rotation curve ^/Bq is raised. The 
result is a rigorous degeneracy of the stellar vs. halo mass 
distribution, insensitive to Hq and the observed time delay 

Atobs- 

3.2. Break the degeneracy from observable flux ratio and 
lens light profile? 

The flux ratio of the saddle image and the minima im- 
age can in princinple differentiate some of the models, as 
shown by Fig. 3. But assuming a reasonable ±0.1 mag 
error with the magnitude, and a 10% error with the effec- 
tive radius, most of the models are in fact degenerate. In 
particular, it will be difficult to constrain the mass of the 
stellar component mo. 

Observations of a well-resolved stellar lens can fix the 
effective radius a = R^/Re and perhaps the cuspiness 
2 — a. This reduces the available lens models drasticly. 
Nonetheless, the degeneracy between lensIII and lensIV 
"iWplies that it is still problematic to differentiate between 
models with dark halo and models without. To break the 
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degeneracy one needs at least an accurate measurement of 
the external shear 71 to 10% level (cf. the parameters of 
lensIII and lensIV in Table 1), perhaps by a combination 
of strong lensing and weak lensing data. 

4. CONCLUSION AND COMPARISON WITH EARLIER 
MODELS 

It is possible to construct many very different models 
with positive, smooth and monotonic surface densities to 
fit the image positions. There are also no extra images. 
These models fit the same images, time delay, Hq and cos- 
mology. Some fit the same lens light profile and image 
flux ratio as well. Hence the models are virtually indistin- 
guishable from lensing data. There are severe degenera- 
cies in inverting the data of a perfect Einstein cross to the 
lens models, even if given the Hubble constant and cos- 
mology. These rigourous findings with analytical models 
arc also consistent with earlier numerical models of Saha 
& Williams (2001) and semi-analytical models of Zhao & 
Pronk (2001). 

Among the acceptable models the rotation curve can be 
Keplerian or flat (Fig. 2a), so lensing data plus Hq cannot 
uniquely specify the lens mass profile. Among models in 
the literature, isothermal models and other simple smooth 
models of dark matter halos of gravitational lenses often 



predict a dimensionless time delay HoAtobs much too small 
(e.g., Schechter et al. 1997) to be comfortable with the ob- 
served time delays Atobs and the widely accepted value of 
Hq ~ 70kms~"'^Mpc~"'^. Naively speaking the high Hq sug- 
gests a strangely small halo as compact as the stellar light 
distribution (Kochanek 2002). But our analytical models 
suggest that there are still many other options. The high 
Ho implies that k is small ~ 0.1 — 0.2 at the images, but 
this does not necessarily imply a rapidly falling density. 
A high Hq does not necessarily mean no dark halo, and 
models with a flat rotation curve does not always yield a 
small Ho (e.g., compare lensHI and lensIV in Figure 2b, 
both satisfy ho = 0.7). We also comment that it will be 
difficult to determine the cosmology from strong lensing 
data alone because the non-uniqueness in the lens mod- 
els implies that the combined parameter w{Ho, i^, zi,Zs) is 
poorly constrained by the lensing data, even if Hq, zi and 
Zs are given. 
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Table 1 

Plausible lens parameters to fit the images of a perfect Einstein cross and time delays with 
Hq = lOO/io = 70kms"^Mpc"^ (cF. equations (15) and (18)). 



Model 


Stars 

Mo 

"r^ 


Half-mass 


Cusp 

2 - a 


Halo 

bo 
Re 


S 


Shear 

71 


Conv. 


KhiRE) 


Mass 

M,{<Re) 
R% 


Mh(<RB) 


lensi 


.2002 


.5 


1. 


.6220 


.2 


-.4445 


.0222 


.1110 


.1334 


.4220 


lensll 


.7110 


1. 


1. 


.2988 


.1 


-.4458 


.0888 


.0494 


.3555 


.1988 


lensIII 


.4901 


.5 


0. 


.2759 


.1 


-.4317 


.0784 


.0379 


.3921 


.1759 


lensIV 


.7129 


.5 


0. 


.0013 


0. 


-.4291 


.1140 


.0006 


.5703 


.0013 
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Fig. 1. — Panel (a): Time delay contours in the XY quadrant in intervals of 3 days for the four lens models. All models reproduce the 
same minima and saddle image positions and time delay. Panel (b) : Cuts of the time delay surfaces -i ^ ^ function of the 

V^saddlo ^mimimai 

radius -R along the radial direction from the lens to the time delay minima (marked by "m" on solid curves to the right) and from the lens 
to the saddle image (marked by "s" on dashed curves to the negative R side) for all four lens models. All models use ho =0.7 and the same 
time delay between the images "s" and "m" . 
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Fig. 2. — Panel (a): The deflection strength of the lens b = M{R)/R as a function of distance from the lens; this is effectively a rotation 
curve of the lens. From top right to down right, from thinner to thicker lines, the models are color coded as red (lensi), blue (lensll), green 
(lensIII) and black (lensIV). Both Keplerian and flat rotation curve models are found. Panel (b): Log-Log plot of the surface density profiles 
k{R) of the four lens models. All models have very low k at the images near —0.3 < logij < 0, consistent with a high Hq. In all models 
except for lensIV (black), stars are dominated by dark matter halos just beyond the Einstein radius at logi? ~ 0. 
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Fig. 3. — The relative magnitude of images as a function of stellar mass of the model mo- From top to down (color coded as pink, cyan 
and orange) the painted polygons correspond to stellar cusps of slope 2 — a = 1, 0.5, 0. Clockwise from the rightmost corner the four corners 
of each polygon correspond to [S,a] = [0,0.55], [0.2,0.55], [0.2,0.5], [0,0.5]. If these models are extended to 5 = 0.3, then all models merger to 
a point on the vertical axis, i.e., a pure-halo model with mo ~ 0. 



